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ABSTRACT: We investigate how Liischer’s relation between the finite-volume energy 
of two pious at rest and pion scattering lengths has to be modified in quenched QCD. We 
find that this relation changes drastically, and in particular, that “enhanced finite-volume 
corrections” of order L® = 1 and i -2 occur at one loop (L is the linear size of the box), 
due to the special properties of the j/ in the quenched approximation. We define quenched 
pion scattering lengths, and show that they are linearly divergent in the chiral limit. We 
estimate the size of these various effects in some numerical examples, and find that they can 
be substantial. 
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1. Introduction 


Almost all lattice QC'D efforts so far have been devoted to the computation of the hadron 
spectrum and hadronic weak matrix elements. However, much more is known about hadrons 
than just their spectrum and weak decays, and one would like to derive this rich phenomenol¬ 
ogy quantitatively from QC'D. The example which concerns us here is the computation of 
pion scattering lengths. 

It is difficult to obtain scattering amplitudes directly in lattice QC'D, since this would 
involve an analytic continuation to the physical regime of numerically obtained Euclidean 
four-point correlation functions. However. Liischer has shown [1,2] that the energy eigen¬ 
values of states with the quantum numbers of two-particle states in finite volume admit a 
systematic expansion in 1/E. where L is the linear dimension of a spatial volume with peri¬ 
odic boundary conditions. In addition, he showed that the coefficients of the powers of 1/7 
in this expansion are related to the infinite-volume elastic scattering phase shifts at values of 
spatial momenta that occur in the finite volume, for energies below the inelastic threshold. 
In particular, for the shift A E in the lowest two-particle energy eigenvalue E for spinless 
bosons with mass m which are at rest, one has 

A E = E - 2m = + ci^ j + 0(L~ 6 ). (1) 

where a q is the l = 0 scattering length, c\ = —2.837297, and c 2 = 6.375183. It is important 
to note that equations such as these are only expected to be applicable in the domain where 
the finite-volume corrections to single particle energies, which are exponentially suppressed 
[3]| are indeed small compared to the right hand side of eq. (1). 

Relations such as these simplify the numerical effort needed to obtain information on 
phase shifts enormously. For instance, one can obtain the / = 0. I = 0 and 1 = 2 pion 
scattering lengths from eq. (1) by extracting the lowest energy levels in the 7 = 0,2 channels 
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from the large time behavior of the Euclidean correlation functions 


Cl= o(f) = |c+_(e - ic++(o, 

0=2 (*) = C++(t), 

C+-(t) = (0|7r + (f)7r-(t)7r-(0)7r + (0)|0) coni 

C++{t) = i(0|7r + (f)7r + (f)7r-(0)7r-(0)|0) co „, (2) 

where 7r ± (f) = fd 3 x 7 r ± (x,f) are zero spatial momentum charged pion fields at time I. 
The subscript “con” means that we should exclude graphs where the fields at time t are 
disconnected from the fields at time 0. The extra factor one half in the last equation comes 
from the fact that the initial and final 7r + 7r + states both need a factor l/\/2 to be properly 
normalized because they art- states with two identical particles. This method has actually 
been used in attempts to compute the 1 = 2 scattering length [4,5], and. more recently, the 
1 = 0 scattering length (as well as the 1 = 2. the pion-nucleon and the nucleon-nucleon 
scattering lengths) [6], 

However, all these computations were done in the quenched approximation, and it is not 
clear a priori that Liischcr's analysis carries over to the quenched approximation without 
modification. New infrared divergences occur in quenched QC'D in the chiral limit [7.8.9], 
and modification of the L dependence of equations like ecp (1) might occur, at least for those 
particles sensitive to the chiral limit. For instance, it was shown that the imaginary part of 
the quenched tt —► tt"^" tt~ scattering amplitude diverges at threshold ( i.e at zero relative 
momentum) when calculated using quenched chiral perturbation theory (ChPT) [10]. (For 
quenched ChPT. see refs. [11,7,9.12], In this paper, we investigate in more detail the case of 
pion-pion scattering at threshold in quenched chiral perturbation theory. We will consider 
only the case of degenerate quark masses. 
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2. Calculation in Quenched ChPT 


The correlation functions C.\ _ (t) and C++(t) can be calculated in the ChPT expansion. 

Formally, with : (g| t) the operator that creates a two-particle state with definite quantum 
numbers, the Euclidean correlation functions C = Cj = q $ can be expressed as 
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+ 0(A 3 ), 
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a . . a a. 

where A is an expansion parameter, proportional to ?n"/(167r /") in unquenched ChPT (we 
will come to quenched ChPT later). We have used the expansions 

E a = 4 0) + a4 1} + a 2 4 2j + 0(A 3 ). 

4L = 4°) + a 4 1} + a 2 4 2) + o(a 3 ). 

The sum over a is a sum over all intermediate states except for the vacuum |0): these states 
are eigenstates of the complete theory (fie., to all orders in A). From eq. (3) we see that the 
perturbative corrections to the energy eigenvalues can be extracted from the terms linear in 
/. after taking out the factors Z^exp(-E^t). We can calculate the energy eigenvalues in 
a finite box A 3 with periodic boundary conditions by restricting all spatial momenta p to 
p = 2x11 /L with n £ 77?'. Calculating to one loop order (fie., to order A 2 ) and expanding 
the result in l/L, we should recover the first two terms of eq. (1). 

The calculation of the 1 = 0 and 1 = 2 correlation functions at tree level in ChPT is 
straight forward, giving 
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from which we obtain the well known result 


n I =0 
r 'o 


777 ; 


ICtt f n 


9 • 


4= 2 


777, 


8- fl 


( 6 ) 


(/tt = 132 Mi 1' in our conventions), using eqs. (1,3). This tree-level result is also valid in 
the quenched approximation. If we extended the calculation to one loop, we would find the 
first two terms of eq. (1). with the infinite-volume one-loop expression for oq in the L ~3 
term, and the tree-level expression for a 0 in the L~ 4 term. For the L~' 3 term one would 
need to go to higher order in the chiral expansion. 


We now wish to determine the quenched approximation version of eq. (1) for 7 = 0 and 
7 = 2 pion scattering. We will show that eq. (1) changes drastically. In particular, "‘enhanced 
finite-volume corrections” to the infinite-volume result E = '2m n occur at order LP = 1 and 
L~-. We will first present the calculation, and then discuss this remarkable result. 


In order to calculate the quenched one-loop corrections to eq. (5), we employ (Euclidean) 
quenched ChPT, which was systematically developed in refs. [7,13]. It was shown there that 
the ?/ meson in the quenched approximation has both single and double pole terms in its 
propagator: 


(?/($>/( 4 )) = Hp + q) 
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where p is the parameter equivalent to the singlet part of the ?/ mass in unquenched QCD. 
From experiment, /r/3 = (500 \h\ - in full QCD: quenched numerical computations get 
roughly this value [14.15.16,17]. The vertex /r can actually have momentum dependence 
[7], which we will ignore for the purpose of this paper. It is the double pole in eq. (7) that 
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will lead to the enhanced finite-volume corrections when an rf appears as an internal line in 
a one-loop diagram contributing to pion scattering. 



(a) (b) 



(c) (d) 


Figure 1. p? -dependent s-channel (a and b) and, i -channel (c and d) contributions to pion- 
pion scattering. Incoming (outgoing) particles are on the left (right). An internal line with 
(without) a cross refers to the double (single) pole term in eg. (7). 

There are two types of contributions containing double pole terms. Figs. 1(a) and (b| are 

the .s-channel meson diagrams for the k~ —» tttt amplitude ( C-\ _) in quenched CliPT, 

and figs. 1(c) and (d) are the f-ehannel diagrams for the same amplitude. C++ gets ^-channel 
contributions from diagrams like figs. 1(c) and (d) as well as crossed (//•channel: versions. 
The crosses denote the /r-vertex corresponding to the second term in eq. (7), see ref. [7]. 
Note that any unusual behavior of quenched correlation functions due to the double pole in 
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eq. (7) is always accompanied by /r dependence. This implies that in order to uncover this 
behavior, we need to calculate the /i" dependence of correlation functions. We will do this 
for C -1 _( t ) and C++(t) to one loop. 

These diagrams have been calculated before in Minkowski momentum space (in infinite 
volume) [10]. and it was found that they diverge on shell for p —> 0 as 1/jfi (with p the 
relative momentum in the center of mass frame): whereas they are perfectly well defined 
in the Euclidean regime. This seems to indicate that no Hamiltonian formalism can be 
developed for the quenched approximation and hence that equations like eq. (3) cannot be 
applied. However we can still parametrize C(t) in Euclidean space as in ecp (3) up to terms 
linear in t, and extract the quantity A E a = A+ A+ .. . from the term linear in 
t inside the square brackets. (The terms quadratic in t, however, will not be of the form 
indicated in eq. (3) [5]). We will take 4-A/A ; as a definition of the two-particle energy in 
the quenched approximation. Note that this prescription coincides with the definition used 
in numerical work [4,5.6]. 


The calculation of the diagrams in fig. 1 is straightforward. Only the mass term in the 
chiral Lagrangian contributes to the 4-meson vertices, since in the degenerate mass case the 

kinetic energy term gives no couplings to the ?/. The .s-channel contribution to C-\ _ (t) 

(tigs. 1(a) + 1(b)) is 
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where o;(k) = y/'m| + k- and k = 2/rn/T with n G 2$. t\ and to are the Euclidean time 
coordinates of the two vertices in tig. 1, and k is the spatial loop momentum. The integrals 
over t\ and to can be carried out, and collecting only the terms proportional to t exp(—2m 7r t) 
we obtain 
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There are also terms proportional to exp( — 2w(k)t| (multiplied by powers of t) arising from 
eq. (8), corresponding (for k / 0) to states with nonzero relative 3-momentum and hence 
higher energy in eq. (3). For k = 0 such terms combine with explicit exp(— 2m 7r f| terms 
to yield contributions finite at k = 0. which have been explicitly separated out in ecp (9). 
From eq. (9) and from figs. 1(c) and (d) and the corresponding diagrams for C++(t), using 
eqs. (2,3), correcting for the renormalization factors = L®/(4?« 2 ) (cf. ecp (5)). and 
noting that all momentum sums are ultraviolet finite, we finally obtain (with A/•.',* , ,n ' loor> = 
>?Ea\ cf. eq. (4)) 


•7-0 l<H,|> = m x EQ(m K L) A + . l|j( ///-/.) be + O 
p.one l°°p _ B-2(m K L) b 2 + Ayfy&irL) be + O 
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where the 1/w 5 and 1/u? terms come from the t- (and u-) channel diagrams. Note that 8 and 
e have to be taken as independent expansion parameters in quenched C'hPT [7,12], and that 
only 8- ( he., /r- ) dependent corrections have been included. If the momentum dependence 
of /.r had been included (/r —» /.r + otp“) then there would be additional terms of order o 6e 

. . A 

in eq. (10). There are ^-independent contributions as well, but those are proportional to f~. 
We will not calculate such terms in this paper. (The 0(e~) terms include contributions from 
the $o-depenclent potentials — sec ref. [7].) There are no Independent contributions from the 
one-loop diagrams with a single six-point vertex, because such contributions are absorbed 
into the renormalization of the pion mass in the tree-level result. This is an example of a 
phenomenon which occurs quite generally for degenerate quark masses [9], 

The complete results to one loop for the energy shifts are then 


Ei =o - 2m* = A E I=0 = APg + A£^ loop , 

E I=2 - 2m* = A£j = 2 = A+ AT"J 2 loop . (13) 


with A E onc loop given by eq. (10) and 


A ptree _ 

AA/_o - 


A = 
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4 /|i 3 
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(14) 


We note that m n in eq. (13) and all subsequent equations is the renormalized, finite-volume 
pion mass, including all one-loop corrections [7]. 


Some momentum sums that appear in ecp (TO) are more singular near k = 0 than in the 
unquenched case. These momentum sums can expanded in L ~ 1 using a result established 
by Liischer in ref. [1]: 
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with integer q > 0 and 0) = /(/.’“)|^. 2 _q- valid up to corrections vanishing faster 

than any power of Z _1 if all derivatives of / are square integrable. A/,. is the Laplacian with 
respect to k and z(q) = Zoo(^.O) is a zeta function, with 

* (3) = 8.40192397. z{2) = 16.53231596. ~(1) = -8.91363292, r(0) = -1. (16) 


The t- and //-channel contributions do not contain any singular momentum sums, and instead 
of ecp (15) we have 
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for the nonsingular case, again up to corrections vanishing faster than any power of Z _1 
Using eqs. (15.17). we obtain large m K L approximations for the coefficients Bq o and Aq^' 
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The Z _0 terms in ecp (18) come from the integrals in eqs. (15.17), and the Lr = 1 and L ~" 
terms in the 7 = 0 channel are the enhanced finite-volume corrections that we announced 
earlier. We also see that the Z - ^ term does not follow the pattern of ecp (1). as can be seen 
from the fact that it is S dependent, whereas the tree-level scattering length is not. There 
are no enhanced terms in the 1 = 2 channel. Note that the corrections to ecp (18) vanish 
faster than any power of Z -1 . 


The enhanced finite-volume “Corrections to A E in the quenched approximation make 
it impossible to define a scattering length via ecp (1). The only alternative is to drop the 
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enhanced terms and define the scattering length to be simply —m n /( 47t) times the coefficient 
of the L ~** term in eq. (1). Using the i ref level result for 0 q ~and eq. (11), we then obtain 
from ecp (10) 


n I=Q - 
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«0 = 


i m 
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8k v 3m. T 

Note that even after removing the enhanced finite-volume corrections, these results still 
diverge in the chiral limit. We have uncovered yet another example of the bad chiral behavior 
of the quenched theory [7.8.9,12], 


For some values of m n L we have computed the values obtained from the definition of 
the coefficients AQ^( m KL). -Bo.2( m kL) (eq. (12)) and also from the approximate expressions 
given in ecp (18): see tables 1 and 2. From these tables, it is clear that, for the smaller 
values of m n L currently used in numerical computations, the exact expressions will have to 
be used. 


3. Numerical Examples 

As a “real world” numerical example, we take m K = 140 MeV, f n = 132 MeV, S = 0.1, 
and m w L = 6. We get 

AE'jZq = -1.3 MeV| AfiJ” loop = -0.3 MeV, 

= 0.36 MeV; A-E™° 2 loop = -0.07 MeV. (20) 

Note that quenched ChPT appears to be working reasonably well, with ~25% one-loop 
corrections to tree-level values. However, because of the enhanced finite-volume corrections, 
quenched ChPT breaks down for larger L. For example, for m n L = 8 we find a 63% one-loop 
correction in the 7 = 0 channel; for m^L = 12, a 240%. correction. In addition, as L gets 
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smaller the sums in eq. (12) get large contributions from small k and are much larger than 
their (already large) asympotic expansions. This means that, as one might expect, quenched 
CliPT also breaks down for small volumes. For example, while m n L = 4 still gives acceptable 
one-loop corrections, for m w L = 2.5 we find a 94% correction in the 7 = 0 channel. 

We note also that 5 as large as ~0.18 is not excluded; 8 = 0.18 is what one would get 
by putting in the full QC'D value for /r: /r/3 = (500 MeV ’. Indeed it is possible that the 
difference between 8 = 0.18 and the 8 ~ 0.1 found in simulations [14.15,16,17] is simply due 
to the larger values of f K which accompany the larger quark masses used on the lattice. With 
8 = 0.18. (and m^L = 6). we have AEj 1 ^ 00 ^ = —1.3 MeV and AEj 1 ^, loop = —0.3 MeV, 
again indicating a breakdown of quenched CliPT. 

For the results given in ecp (10) to be useful in numerical simulations, the terms linear 
in t in eq. (3) must not be obscured by terms with higher powers in t. Furthermore, one 
must be able to separate numerically the contribution of the lowest two-pion state (with 
relative momentum k = 0) from that of the closest excited states (with relative momentum 
|k| = 2n/L). We therefore demand 

|A£ troc #| < 1, (21) 

and 



It is not hard to satisfy these conditions in our “real world” example with m n L = 6. Putting 
|A/-.'" V< ' t\ = 0.1 we find Qj_q = 10 and Q/_ 2 = 35. Note that, because A7f tree ~ /. :i . 
fl actually decreases for smaller m K L if |A7% reo t\ is held fixed. For example, af m^L = 4, 

^ 7=0 = 3 - 

We can also calculate the size of the one-loop quenched CliPT effects in a recent numerical 
computation of scattering lengths by Kuramashi et al. [6]. They used a 12'^ x 20 lattice at 
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(3 = 5.7, and considered pion masses of 0.29 (using staggered fermions) and 0.508 (using 
Wilson fermions), in lattice units. The values of f K were 0.187 and 0.143 respectively [18] 
(our / vT is \/2 times that of ref. [18]): we take 8 = 0.1. 


In their numerical computation Ivuramashi et al. did not include the “double-annihilation" 
diagram (in the terminology of ref. [5]), which means that in comparing with eq. (10) we 
should drop those 0(8“) terms that come from -4-channel amplitudes. This leads to a modi¬ 
fication of the expression for /ipi /»-/.; ( cf. eqs. (10.12)). which now becomes 


B ( 


207T 4 ?77 
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For the various one-loop Contributions to A F we obtain using the lattice values for f K (in 
lattice units): 



aAE I=Q 

. j-ono loop 

flAjfc T=0 

aAE^l 

» r-iOiie loop 
Cl LA Hi j _ c) 

n in K = 0.29 

-0.029 

0.0002 

0.0083 

-0.0017 

am 7 t = 0.508 

-0.050 

0.005 

0.014 

0.001 


W T e see that in this case the one-loop corrections are quite small — anomalously so in 
some cases because of cancellation between the 8 and he terms. (For S = 0.18, the one-loop 
terms are less than 20% of the tree-level terms, except for AEj 1 ^ loop which is of the same 
size as A Ej r S% for mil- = 0.29.) On the other hand, the conditions ecp (21) and eq. (22) 
are just barely satisfied, if at all. in the range of t (4 to 9) in which they fit. For example, 
for arn^ = 0.29, [AFj^q t\ = 0.26 at 7 = 9, and Q = 2.4 at t = 4, For mil- = 0.508. 
|A/•.[’)[] i | = 0.45 at i = 9, and 0 = 1.8 at t = 4. Since there are six excited states with 
|k| = 2 /t/Z,. considerably larger values of Q would be needed to be confident that the excited 
states are not contaminating the results. 


. O . / ■ 

Removing the 8“ terms in ecp (19) that come from s-channel contributions, we would 
obtain for the 1 = 0 scattering length (cq = “ remains unchanged): 

(24) 
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For am K = 0.29 the f>~ term leads to a 3% correction, and for am % = 0.508 it is 1%. These 
numbers are 11% and 2% respectively in the 1 = 2 case. 

4. Origin of the Enhanced Finite-Volume Corrections 

It is not difficult to understand the origin of the enhanced finite-volume corrections to 
eq. (1) in the quenched theory. First, however, we need to have an intuitive picture in 
a normal theory of how the L dependence of the right hand side of ecp (1) arises within 
Euclidean relativistic perturbation theory. (See also ref. [5].) 

We begin by examining the factors of L and t in the zeroth order contribution to the 
correlation functions in eq. (2). At this order, we just have a disconnected product of two 
separate pion propagators from time 0 to time t. Since each of the four external fields is 
being integrated over space, there would be a factor of L were it not for the fact that the 
propagator falls exponentially for |.ij > 1 /f» vT , which forces the start and end of each line 
to be close to each other in space, and reduces the L 1 - to L°. Furthermore, each pion line 
contributes a factor of e~ m ^; together they produce the factor which appears in 

ecp (5). 

At next order, “tree level,” we have one interaction vertex. This vertex is integrated 
over all space-time., but now all 5 points (4 external plus the vertex) must be spatially close 
to each other and only a L? survives. The time, I 1 , of the vertex can however be anything as 
long as t > t 1 > 0, so this gives a factor of — in other words, this is a contribution 

to the two-particle energy. Relative to the zeroth order term, this contribution is suppressed 
by L ~^ and therefore contributes to the first term on the right hand side of ecp (1). 

On the other hand. if. for example, t 1 > t. the total time for propagation is now t+2(f' —t ). 
and the contribution is suppressed by an additional factor r - E. So only times f 1 > f 

with t' — t ~ 1/m 7 T contribute, and this region of integration over i' produces just a factor 
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~ (l/m 7T )e . in other words this is a “Z factor 1 ' correction in the sense of eq. (3). not 
an energy contribution. 

Now let us go to the one-loop .s-channcl diagram (in the full theory, or any normal theory, 
like (p 4 ). Call the times of the two vertices t\ and t- 2 . with the vertex closest to the outgoing 
( t ) lines. In space, all six points (4 external plus 2 vertices) must be close (~ 1 /m -: to each 
other, so we start with a factor of Z 3 . Consider the integration over t\ and t-j- If t > t\ = 
t -2 > 0, the integrand goes like . However, as t\ moves away from t 2 ■ (but still in the 

order t > t\ > t^) the integrand is suppressed by an additional factor fi _ -Mk)-r»7r)(Z -U ). 
where k is the momentum mode in the loop. This forces t\ — t >2 ~ tj. = l/(w(k) — m^). 
In infinite volume, the lowest scale in the loop would Ire k ~ m w . For such k. tj, ~ 1 / in-. 
t\ to. and the integration over t\ and t .2 produces single factor of t. This is the normal 
one-loop contribution to the energy of order L~ 3 relative to the disconnected diagram. In 
other words this is just a correction to oq in the L~ 3 part of eq. (1). 

However, in finite volume one must sum over a discrete set of k's. and a lower scale for k. 
1/Z. is available. Indeed, to estimate the difference between the sum over k and the integral 
over k. one can simply look at the contribution of the lowest nontrivial momentum modes, 
k ~ 1/Z. (For the contribution of k = 0. sec below). 

First note that the factor Z 3 ( which one would have in “infinite volume” is re¬ 

placed by (>0^ in finite volume, so the contribution of any single finite-volume mode is down 
by L~ 6 from the zeroth order term. The modes k ^ 1/Z have have cj(k) — m vT ~ Z _ “, so 
Tj, ~ Z". Now the integration over t\,t2 from the region / > / | > 1 2 > 0 gives r^j T k tr 2m *% 

ty # _ a 

Since tj. ~ Jr this gets boosted back to Z relative to the zeroth order diagram. In other 
words, the difference between the sum and integral over k produces the OqL~ 4 term in the 
two-particle energy. Contributions where for example t -2 > t or t\ > t. etc., arc suppressed 
with no gain in Z factors, so these give higher order corrections in l/Z. 

When k = 0, there is no suppression at all for t > i\ >t 2 > 0: t\ and t^ move freely in 
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this range. But then the integration over ti-i-i gives a factor f", so this is not a contribution 
to the one-loop energy (it is an iteration of the tree-level energy shift). “Out of order" 
contributions are suppressed by c~^ m ' K ' r , where r / is the amount of time out of order. These 
can give but there is no L enhancement, so these are terms of order T -6 in the 

energy. 

Note that the t- and (/-channels do not contribute L~ 4 terms in the energy: the internal 
lines always contribute e ^^(k)(f] ) ( no ^(k) — m K in the exponent). So tj, ~ 1 jrn^ always, 

and there is no boost in the power of L from the integrations over o. 

The quenched “anomalous terms are now easy to understand. An “x" (f,r vertex) on 
an internal line can be at any time with no additional suppression as long as t\ > t x > t-2- 
So the integration over t x produces an extra factor of 77 . Therefore each “x” gives two more 
powers of L. starting from the term L~ 4 in the energy. The single-x diagram's contribution 
to the energy is thus of order L~“, and the double-x's, order 1 . Again, there are no such 
enhanced contributions in the t- and //-channels. 

We conclude this section with the remark that, in the quenched theory, the enhanced 
terms in 1/L are generated exclusively by the lowest modes in momentum space, i.e.. modes 
with n ~ 1, where k = 2~n/ L. The sum over n for the enhanced terms converges well before 
n ~ w n L. On the other hand, the L~ 4 terms in the energy in a “normal" theory ( e.g ., full 
CliPT). which are analogous to the enhanced quenched terms, are generated by all modes 
up to n ~ m^L, not just by the 11 ~ 1 modes. This is because the enhancement by L~ from 
phase space for modes with n ~ ///. - /. competes with the Tr enhancement from 77 . for modes 
with n ~ 1. Indeed, in a normal theory, the sum over all modes with n</?? vT i contributes 

O A 

at order Z, -J to the energy, and the terms arise from the difference between the sum 
and the integral over such modes. 
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5. Conclusion 


We have shown how to adapt Liischers analysis of the relation between the finite-volume 
energy shifts of two-particle states and scattering lengths to the quenched approximation 
of QCD. These energy shifts are defined directly from the quenched Euclidean correlation 
functions, in a way consistent with the way it was done in recent numerical computations 
of two-pion correlation functions. We have calculated the finite-volume energy shifts to one 
loop for the case of two pious at rest in quenched chiral perturbation theory; our main result 
is contained in eq. (10)-eq. (14), 

The one-loop corrections are very different from those of the full theory. In particular, 
enhanced finite-volume corrections” of order 1 and Jr ~ appear, to be compared with the 
tree-level contributions which are of order Jr ''. This indicates yet another breakdown of 
quenched chiral perturbation theory (and presumably quenched QCD) in the infrared, and. 
as in other quantities, this behavior originates from the special role of the // in the quenched 
approximation. We have examined the size of these effects in some examples, including one 
in which we take the pion mass, the pion decay constant and the volume from the most 
recent numerical computation of the energy shifts [18]. The range of applicability of the 
whole analysis as a function of the volume, the pion mass, and the // mass is also discussed. 

Because of the completely different nature of the quenched one-loop corrections (as 
compared to the unquenched case), and in particular because of their infrared divergent 
behavior, the quenched two-pion energy shifts will only be close to the full QCD results 
if the one-loop corrections are small compared to the tree-level terms. A consequence of 
this is that one can expect agreement of quenched computations with full QCD at best to 
25%, which is the size of one-loop corrections to tree-level chiral perturbation theory in the 
unquenched case. (See for instance a recent compilation of tree-level and one-loop chiral 
perturbation theory results for pion scattering lengths by Gasser [19].) 
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Presumably no satisfactory Hamiltonian formulation exists for quenched QCD. Indica¬ 
tions in this direction are the enhanced finite-volume corrections uncovered in this paper, 
as well as the (related) fact that the pion scattering amplitude, continued to Minkowski 
space, diverges at threshold [10], This implies that the energy shifts and scattering lengths 
as defined in this paper from Euclidean correlation functions do not have a direct physical 
meaning. One may nevertheless hope that the quenched values for these quantities are close 
to similarly defined quantities in full QCD (at least for some range of parameters), which 
do have a physical interpretation. As discussed above, a necessary condition for this is that 
the values of these quantities are given essentially by tree-level chiral perturbation theory. 
Finally, we would like to remark that a nonperturbative analysis along the lines of ref. [2] 
does not seem to be possible, since for that approach a Hamiltonian framework is indispens¬ 
able. The only analytic handle we have on pion physics in the quenched approximation is 
Euclidean quenched chiral perturbation theory. 
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Table 1. Bq and Aq as a function of m n L. The third (fifth) column (-Bo exp/ Bq resp. 
Aq exp/^o) gives the ratio of Bq (Tq) calculated from cq. (18) and eq. (12). 


m^L 

Bq 

Bq exp /Bq 

-'() 

4^0 expA^O 

1.0 

-3592.1 

1.0202 

9363.3 

0.50044 

1.5 

-315.47 

0.96660 

828.75 

0.50779 

2.0 

-56.287 

0.89612 

151.71 

0.53075 

2.3 

-24.443 

0.84633 

67.895 

0.55550 

2.5 

-14.902 

0.81036 

42.560 

0.57707 

2.7 

-9.4693 

0.77256 

27.993 

0.60252 

3.0 

-5.1370 

0.71373 

16.180 

0.64664 

3.3 

-3.0011 

0.65514 

10.180 

0.69511 

3.5 

-2.1797 

0.61837 

7.7905 

0.72799 

3.7 

-1.6310 

0.58529 

6.1387 

0.76007 

4.0 

-1.1149 

0.54625 

4.5037 

0.80478 

4.5 

-6.7984 

0.51968 

2.9774 

0.86622 

5.0 

-0.48446 

0.54330 

2.1613 

0.91025 

5.5 

-0.39019 

0.60041 

1.6717 

0.93986 

6.0 

-0.34224- 

0.66856 

1.3502 

0.95923 

6.5 

-0.31697 

0.73276 

1.1238 

0.97182 

7.0 

-0.30343 

0.78697 

0.95588 

0.98007 

7.5 

-0.29622 

0.83043 

0.82635 

0.98556 

8.0 

-0.29257 

0.86450 

0.72343 

0.98929 

9.0 

-0.29054 

0.91166 

0.57054 

0.99372 

10.0 

-0.29151 

0.94055 

0.46300 

0.99604 

11.0 

-0.29347 

0.95875 

0.38388 

0.99736 

12.0 

-0.29566 

0.97058 

0.32374 

0.99817 
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Table 2. Bo and Ao as a function of m n L. The third (fifth) column (Bo exp / Bo resp. 
A 2 exp/M) gives the ratio of Bo (Ao) calculated from cq. (18) and eq. (12). 


m n L 

Bo 

Bo exp/-^2 

M 

'2 exp M 2 

1.0 

-1948.2 

0.0067546 

4680.4 

0.016870 

1.5 

-171.08 

0.022791 

413.47 

0.056581 

2.0 

-30.496 

0.053939 

75.145 

0.13134 

2.3 

-13.219 

0.081819 

33.310 

0.19482 

2.5 

-8.0397 

0.10476 

20.682 

0.24433 

2.7 

-5.0891 

0.13137 

13.432 

0.29864 

3.0 

-2.7327 

0.17835 

7.5685 

0.38638 

3.3 

-1.5677 

0.23358 

4.6049 

0.47712 

3.5 

-1.1178 

0.27459 

3.4316 

0.53665 

3.7 

-0.81582 

0.31845 

2.6253 

0.59376 

4.0 

-0.52955 

0.38828 

1.8345 

0.67251 

4.5 

-0.28357 

0.50926 

1.1104 

0.78032 

5.0 

-0.16838 

0.62524 

0.73624 

0.85795 

5.5 

-0.10892 

0.72620 

0.52135 

0.91027 

6.0 

-0.075482 

0.80712 

0.38719 

0.94408 

6.5 

-0.055213 

0.86787 

0.29781 

0.96541 

7.0 

-0.042100 

0.91130 

0.23521 

0.97868 

7.5 

-0.033138 

0.94130 

0.18964 

0.98688 

8.0 

-0.026731 

0.96153 

0.15547 

0.99193 

9.0 

-0.018349 

0.98377 

0.10864 

0.99694 

10.0 

-0.013249 

0.99325 

0.079048 

0.99884 

11.0 

-0.0099145 

0.99722 

0.059348 

0.99956 

12.0 

-0.0076241 

0.99886 

0.045700 

0.99983 
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